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Abstract 

The self-energy and vertex QED radiative corrections (~ Za 2 f(Za)) are 
shown to give a large negative contribution to the parity nonconserving (PNC) 
amplitude in heavy atoms. The correction —0.7(2)% found for the 6s — 7s 
PNC amplitude in 133 Cs brings the experimental result for this transition into 
agreement with the standard model. The calculations are based on a new 
identity that expresses the radiative corrections to the PNC matrix element 
via corrections to the energy shifts induced by the finite nuclear size. 
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It has been discovered recently that there exists a consistent deviation of experimental 
data on parity nonconservation (PNC) in atoms from predictions of the standard model. 
This paper demonstrates that this contradiction is removed by the self-energy and vertex 
QED radiative corrections, which prove to be much larger than anticipated. The corrections 
are evaluated with the help of a new identity that expresses them via similar radiative 
corrections to the energy shifts induced by the finite nuclear size. 

Experimental investigation of the 6s — 7s PNC amplitude in 133 Cs initiated by Bouchiat 
and Bouchiat [p]], was carried further by Gilbert and Wieman 0, and by Wood et al 
who reduced the error to 0.3%, sparking an interest in the atomic PNC calculations that 
are crucial for the analysis of the experimental data. Accurate previous calculations of Refs. 
P]j20[] have recently been revisited by Kozlov et al || and Dzuba et al @J|. Bennett and 
Wieman || analyzed the theoretical data comparing it with available experimental 

data on dipole amplitudes, polarizabilities and hyperfine constants for Cs, and suggested 
that the theoretical error for the PNC amplitude should be reduced from 1% to 0.4%. It 
has been recognized recently that several, previously neglected phenomena contribute at the 
required level of accuracy. Derevianko ]nj found that the Breit corrections give —0.6%, the 
value confirmed in |IT|,|6]|. Sushkov ]n| pointed out that the radiative corrections may be 



comparable with the Breit corrections. The calculations of Johnson et al [I3| demonstrated 
that indeed, the QED vacuum polarization gives 0.4%, the value confirmed in |^|7|,|T5|] . 

Ref. |J indicates that there is a 2.3a deviation of the weak charge Qw extracted from the 
atomic PNC amplitude || from predictions of the standard model ||16|| . More recent works 



i~3|,|8|1, in which the Breit corrections (—0.6%) and the QED vacuum polarization (0.4%) were 
included, give similar deviations 2.2a and 2.0<r respectively. We show that this contradiction 
is removed by the self-energy and vertex radiative corrections. The corrections of this type 
were considered previously by Marciano and Sirlin jT7 and Lynn and Sandars [[TJ| using the 
plane wave approximation that resulted in a small value ~ 0.1%. The expectation of Ref. 
|T3[ is that the Coulomb field of the atomic nucleus should not produce any drastic effect 
on these radiative corrections, which should remain small. This assessment is supported 



by Ref. [14] that mentions in passing preliminary results of calculations indicating that the 



self-energy and vertex corrections altogether are small. The opposite conclusion of Ref. [|7| , 
namely that the self-energy correction may give a substantial contribution, was not decisive, 
because the model approach pursued in this work was not gauge invariant. 

Let us show that there is a precise relation that expresses the QED radiative corrections 
to the PNC matrix element via similar radiative corrections to the energy shifts of the atomic 
electron induced by the finite nuclear size (FNS). This relation can be presented as 

<5pNC,sp = - (#FNS,s + ^FNS,p) , (1) 

where <5pnc,sp is the relative radiative correction to the PNC matrix element between s\/2 
and pi/2 orbitals. The term relative correction used above indicates that the correction is 
divided by the matrix element itself 

, (yy/ 2 |ffpNC# P ,i/2) rad . . 

CpNC,s P — y- [77 r; r — , \£) 

Similarly <5fns,s an d <5fns,p are the relative radiative corrections to the FNS energy shifts 
£tns,s; -Efns,p; for the the chosen si/ 2 and pi/ 2 electron states 
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<5fNS,Z — ^PNS,j/-^FNS,2 ) I — s iV 



(3) 



The operator H PNC in @ describes the PNC part of the electron Hamiltonian induced by 
the Z-boson exchange 

#p N c = (2v / 2)- 1 G F g w p{r) 7s • (4) 

Here G P and Qw are the Fermi constant and the nuclear weak charge, and p(r) is the 
nuclear density. The FNS energy shifts can be presented as matrix elements of the potential 
<5VFNs( r ), which describes the deviation of the nuclear potential from the pure Coulomb 
one, -Efns,; = (^u^I^^fnsI^U^) >^ = s iP- Equality (]l|) may be established for the sum of 
all QED radiative corrections (~ Za 2 f(Za)), or specified for any gauge invariant class of 
them. We concentrate our attention on the self-energy and vertex corrections in the lowest 
order of the perturbation theory described by the Feynman diagrams in Fig. |T|, calling them 
the e-line corrections, though the vacuum polarization is also briefly discussed below. The 
vertex in diagrams of Fig. ([I]) originates from the PNC Hamiltonian (^), the left and right 
external legs describe the wave functions ip s ,i/2( r ) an d ip P) i/2 ( r ) for the considered atomic 
states. The sum of these diagrams is gauge invariant, though each one of them is not. We 
can use this fact for our advantage, choosing a gauge in which the vertex correction Fig. 
|l](a) is zero. To see that this is possible we can consider, for example, the gauge in which 
the photon propagator is D^ y {k) = {g^ v — S,f{k)k^k u )/k 2 , where f(k) is some nonzero, 
nonsingular function of k, and tune one parameter £ to annihilate the vertex correction for 
the chosen pair of atomic Si/2,Pi/2 levels. In this gauge only the self-energy corrections Fig. 
|l|(b,c) contribute to the PNC transition between this pair of states. 

The PNC Hamiltonian (|J) is localized inside the nuclear interior r < R, where R is the 
nuclear radius. In contrast, the radiative corrections take place at separations comparable 
with the Compton radius r ~ mr x which is much bigger than the nucleus, mR <C 1 (rela- 
tivistic units h = c = 1 are used, if not stated otherwise). This difference of the two scales 
allows us to simplify the problem. Consider the diagram Fig. |l|(b) in which the initial wave 
function (the left leg) is ^,1/2 ( r ) • The radiative correction induced by the self-energy oper- 
ator results in a variation of this wave function S Iad ^) St i/2(r) which we need to evaluate at 
the nucleus, where the weak interaction takes place. In this region the shape of the function 
does not change, because the perturbation caused by the radiative correction is localized far 
away. Thus, inside the nucleus £ rat V.s,i/2( r ) = C T s ajd ip s ,i/2( r ), where C r s ad is an r-independent 
factor. Similarly, for the right leg (diagram Fig. |l](c)) 5 ta ' d ipp,i/2(j t: ) = Cp ad ip P ,i/2 ( r )- Using 
these variations of the wave function we express the relative radiative correction to the PNC 
matrix element (0) in terms of the factors C^ ad , Cp ad 

(>PNC,sp = <^ s + <^ p ■ (Oj 

We took into account here that in the chosen gauge (^s,i/2|-£fpNc|^p,i/2) rad = 

(5 r ^V/2|#PNc|VW/2> + (^,l/2|^PN C |5 rad ^,l/2)- 

Let us discuss now the e-line radiative corrections to the FNS energy shift for the S1/2 
level that are described by the same Feynman diagrams in Fig.[l| in which the vertex is 
given by deviation of the potential from the pure Coulomb one <5VpNs( r )- We can again 
choose the gauge in which the vertex correction is zero for the given S1/2 state. Moreover, 
we can assume that the vertex radiative corrections are zero simultaneously for the FNS 
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energy shift and for the PNC matrix element (for the chosen pair of si/2,Pi/2 states). This 
is possible because gauge transformations include an infinite number of parameters. The 
gauge D tll/ {k) = (g^ u — f(k)k^k l/ )/k' 2 presents them via an arbitrary function f(k) that can 
be chosen to satisfy the conditions formulated above. In this gauge the radiative correction to 
-Efns.s = (V's,i/2| ( 5Vfns|V ; s,i/2) is expressed via the variation of the wave function 5 rad ^ Sj i/ 2 (r), 
which is essential only inside the nucleus where <5VFNs( r ) is located. Thus, using arguments 
similar to the ones that led us to (|5|), we find the relative correction to the FNS energy shift 

£fns, s = 2C™ d , (6) 

where the factor 2 accounts for two diagrams (b) and (c) in Fig.|l| that give identical contri- 
butions. Similarly for the px/ 2 level, 

(Ws,p = 2Cp ad . (7) 

The factors C^ ad , C™ d are determined by the radiative corrections, being independent on 
the nature of a perturbative operator localized on the nucleus. They have, therefore, the 
same values for the PNC matrix element (^) and FNS energy shifts (§|),(0). Combining these 
equations we immediately derive Eq.(|l|). The only parameter which governs the accuracy 
of the presented derivation is the smallness of the nucleus compared with the Compton 
radius. This makes Eq.(|l]) a very accurate identity. Similarly we derive this identity for 
the contribution of the QED vacuum polarization. The result can be compared with Ref. 



15| that presents explicit variations of S\/2 and p\/2 wave functions at the origin induced by 



the vacuum polarization (see Eq.(43) of Ref. |15|). Using these wave functions to calculate 



corrections to the PNC matrix element and FNS energy shifts, we find perfect agreement 
with Eq.flD. 

Note that we do not consider here radiative corrections of the order ~ ct/ir which appear 
in the plane wave approximation. These contributions have been included into the radiative 
corrections to the weak charge Qw (and the renormalization of the charge and electron mass 
in the case of FNS energy shifts). Correspondingly, we subtract the contribution of the plane 
waves from Eq.(^), considering only the part of the corrections that depends on the atomic 
potential ~ Zce 2 f(Za). For heavy atoms this subtlety is insignificant numerically because 
the considered Z-dependent part of the correction is bigger than the omitted Z-independent 
one, as we will see below. 

Eq. ([!]) presents the e-line corrections to the PNC matrix element, which are difficult 
to calculate, in terms of the corrections to the FNS energy shifts that have been well- 



studied both numerically, by Johnson and Soff [IE], Blundell ||20|| , Cheng et al |HJ and 
Lindgren et al [f22]| , and analytically, by Pachucki [23] and Eides and Grotch |24[]. Ref. 
Plfl presents the e-line radiative corrections to the FNS energy shifts for lsi/ 2 , 2si/ 2 and 
2pi/ 2 levels in hydrogenlike ions with atomic charges Z = 60, 70, 80, 90. Eq. (jl]) contains 
relative corrections, therefore we need to calculate the FNS energy shifts E FNS . We did 
this by solving the Dirac equation with the conventional Fermi-type nuclear distribution 
p(r) = po/{l + exp[(r — a)/c}}. Parameters a,c were taken the same as in pi| , namely 
a = 0.523 fm and c chosen to satisfy i? rms = 0.836A 1 / 3 + 0.570 fm. Using the results of 



[ 2~T| and this calculation we obtained the relative radiative corrections shown in Fig. 



In order to include the interesting case Z = 55 and to account for all values of 55 < 
Z < 90 we used interpolation formulae presented in pi]. The relative corrections for the 
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Is and 2s levels are roughly the same size. This indicates that the radiative processes 
responsible for the correction take place at separations smaller than the K-shell radius, 
r < Zam" 1 , which is consistent with the assumption r ~ m _1 above. For these separations 
we can assume that, firstly, the screening of the nuclear Coulomb field in many-electron 
atoms does not produce any significant effect, and, secondly, the relative corrections does 
not depend on the atomic energy level because for small separations all atomic nsi/2-wave 
functions exhibit similar behaviour. These arguments remain valid for the py2 states as well, 
permitting us to presume that the results shown in Fig. || for the 2si/ 2 -levels and 2p 1 / 2 - 
levels of hydrogenlike ions remain valid for S1/2 and pi/2 states of the valence electron in a 
many-electron atom. We obtain the e-line radiative corrections for the PNC matrix element 
using Eq.(|l]) that expresses them via the found corrections to the FNS energy shifts. The 
found PNC corrections, presented in Fig. |2| by the dotted line, are negative and large (much 
larger than the neglected Z- independent part of the corrections). 

Let us discuss the implications for the 6s — 7s PNC amplitude in 133 Cs. The standard 
model value for the nuclear weak charge for Cs ||16| is 

Qiy( 133 Cs) = -73.09 ± (0.03) . (8) 

Ref. [El refined previous calculations of Ref. M extracting from the experimental PNC am- 
plitude of Ref. the weak charge 

Q w ( 133 Cs) = -72.18 ± (0.29) cxpt ± (0.36) theor , (9) 

with the theoretical error 0.5%. It is consistent with Q w ( l33 Cs) = —72.21 ± (0.28) exp t ± 
(0.34)theor that was adopted in |13| by taking the average of the results of Refs. , and 



accepting the theoretical error 0.4% proposed in ||. The weak charge in Eq.(|) deviates 
from the standard model (||) by 2. Oct. 

The e-line radiative correction derived from results presented in Fig. |^ is — 0.73±(0.20)%, 
the error reflects the uncertainty of the radiative corrections to the FNS energy shift for the 
2pi/2 level in Cs -EpNS,2p,i/2 = —0.0001(1) eV p7j (more accurate value can, probably, be 
obtained by extrapolation of much more reliable higher- Z results shown in Fig. |2|). Eq.(|9]) 
combined with the e-line correction gives 

Q w ( 133 Cs) = -72.71 ± (0.29) cxpt ± (0.39) thcor , (10) 

which brings the PNC experimental amplitude of |3[] within the limits of the standard model 
(§). For heavier atoms the e-line corrections become larger, while the error diminishes. For 
the Tl atom, which presents another interesting for applications case, we find the e-line 
correction —1.61%. 

Relations similar to ([!]) can be derived for any operator which is localized at distances 
smaller than the Compton radius. One can even try to apply it to the case of the hyperfine 
interaction (HFI), which has been thoroughly investigated previously, see e.g. [^,^] and 
references therein, though the HFI has a long-range tail ~ 1/r 3 that presents an obstacle 
for our method. However, if convergence of the HFI matrix elements is fast, the relation 
<5fns,s ~ ^hfi s should hold. Here 5^ Fl s is the radiative correction to the HFI for s-levels, 
the primed notation indicates that the Z-independent Schwinger term a/(27r) should be 
excluded (for heavy atoms this subtlety is not important.) Fig.0 shows the e-line contribution 
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to 5ft FI s that was extracted from [[25] using interpolation for all considered values of Z. It 
agrees semi-quantitatively with #fns,s> deviation is less than 33 %. Overall, we observe 
that 5fns,is, <^fns,2s, ^fns,2 P , and 5hfi a a ^ exhibit similar behaviour, they are negative and 
large regardless of the perturbative operator considered and quantum numbers of the wave 
functions involved, which is in line with the main arguments of this paper. 

In conclusion, large QED self-energy and vertex corrections to the parity nonconservation 
amplitude in heavy atoms reconcile the experimental results of Wood et al |5| in Cs with 
the standard model. 

This work was supported by the Australian Research Council. The authors are thankful 
to K.T.Cheng for his calculations of the self-energy corrections to the FNS energy shifts in 
Cs Ref. |3. 
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FIGURES 



FIG. 1. The QED vertex (a) and self-energy (b),(c) corrections to the PNC matrix element 
and the FNS energy shifts. For the electron- nucleus PNC interaction the vertex is given in Eq. (Q), 
for the FNS correction the vertex is produced by the short-range potential that describes the 
spreading of the nuclear charge. 

FIG. 2. The relative radiative corrections (in %) induced by the diagramms of Fig.l. Cor- 
rections to the FNS energy shifts for lsi/ 2 , 2s]/ 2 , and 2p 1 / 2 levels extracted from [21] as discussed 
in the text are shown by thick, thin, and dashed lines. Dotted line - predictions of Eq.(||) for the 
radiative corrections to the PNC matrix element. Dashed-dottet line - relative correction to the 
HFI taken from [25]. 
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Fig.l 




Fig.2 



